Supergranules as Probes of Solar Convection Zone Dynamics 
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ABSTRACT 

Supergranules are convection cells seen at the Sun's surface as a space filling pattern of 
horizontal flows. While typical supergranules have diameters of about 35 Mm, they exhibit a 
broad spectrum of sizes from ~ 10 Mm to ~ 100 Mm. Here we show that supergranules of 
different sizes can be used to probe the rotation rate in the Sun's outer convection zone. We find 
that the equatorial rotation rate as a function of depth as measured by global helioseismology 
matches the equatorial rotation as a function of wavelength for the supergranules. This suggests 
that supergranules are advected by flows at depths equal to their wavelengths and thus can be 
used to probe flows at those depths. The supergranule rotation profiles show that the surface 
shear layer, through which the rotation rate increases inward, extends to depths of ~ 50 Mm 
and to latitudes of at least 70°. Typical supergranules are well observed at high latitudes and 
have a range of sizes that extend to greater depths than those typically available for measuring 
subsurface flows with local helioseismology. These characteristics indicate that probing the solar 
convection zone dynamics with supergranules can complement the results of helioseismology. 



Subject headings: Sun: convection, Sun: rotation 



INTRODUCTION 



Supergran ules were discovered in the 1950s by 



(^ 5 — 10 Mm) have interm ediate lifetimes of 



2 hr ( November et al. Ill98lh . 



Dupergran uies were aiscoverea m tne i9oUs Dy th, , , ■ , r j.i, i j.j. 

TT f j ^^Z 1 1 1 1 . , . , Ihe rotation rate o i the supergran ule pattern 
Hart I ( 19541) but were best cha racterized m the ^ . j > Fn! iTTJinonK \. 



1960s by iLeighton et al. (11962^ w ho gave them 



their name. ILeighton et al. I ( 1962 ) showed that 



this cellular pattern of horizontal flows covers the 
solar surface and that the boundaries of the cells 
coincide with the chromospheric/magnetic net- 
work. Typical supergranules have diameters of 
~ 35 Mm and maximum flow speeds of '^ 500 
m s~^. While the kinetic energy spectrum has 
a distinct peak at wavelengths of ^ 35 Mm, the 
spectrum includes cells at least three times larger 
and extends to much smaller cells where the su- 
pergranule S2ectrimi__bleiids into the granulation 



spectrum (JHathawav et al. 1120001) 



Larger cells live longer than smaller cells. 
Typical supergranules with diameters of ^ 30 
Mm live for ^ 24 hr (jSimon fc Leighton I IT963 
IWang fc Zirin I 119891 ). Typical granules with 
diameters of ~ 1 Mm only live for ~ 5 min 



correlated the Doppler velocity pattern from equa- 
torial spectral scans obtained over several days. 
He found that the pattern rotates about 3% faster 
than the photospheric plasma and faster rates are 
found for the 24- hr time lags from day-to-day than 
for the 8-hr time lags from the beginning to end of 
an observing day. He concluded that larger cells 
dominate the longer time lags and that the obser- 
vations are consistent with supergranules embed- 
ded in a surface shear layer in which the rotation 
rate increases with depth. 

The pre sence of this shear layer was first sug- 
gested by iFoukal fc Jokipii I (|l975[ ) as a conse- 
quence of the conservation of angular momen- 
tum by convective elements moving inward and 
outward in the near sur face layers. Global he- 



( Title et al. 1 119891) . Cells of intermediate size 



lioseismology inver sions ( Thompson et al. I Il996t 



Schou et al. Ill998l) indicate that the shear layer 



extends to depths of 35-50 Mm near the equator 
but may disappear or reverse at lat itudes above 
55°. While lo cal helioseismolog y (JBasu et al 



19991: ICorbard fc Thompson 120021) does not probe 
as deeply, it produces similar results which suggest 
that the rotation rate may not continue to increase 
inward at higher latitudes. 

(l2nnnl i 



Beck fc Schou 



measured the rotation 
of the supergranule pattern using a Fourier tech- 
nique with space-based Doppler data from the 
ESA/NASA Solar and Heliospheric Observa- 
tory (SOHO) Michels on Doppler Imager (MDI) 



(jScherrer et al. 1 119951 ). They mapped the data 



onto heliographic coordinates, took the Fourier 
transform in longitude of data from equatorial 
latitudes, and then the Fourier transform in time 
of those spectral coefficients over 6 10-day inter- 
vals in 1996 which had continuous coverage at a 
15-min cadence. They found that the larger cells 
do indeed rotate more rapidly than the smaller 
cells, but with rotation rates that exceeded the 
peak internal rotation rate at the base (50 Mm 
depth) of the surface shear layer as d etermined 
from global helioseismology (jSchou et al. 1998). 
This discovery led them to conclude that super- 
granules must have wave-like properties in order 
to rotate faster than the flows they are embedded 



However. iHathawav et al. I (20061) showed that 
line-of-sight projection effects on a rigidly rotat- 
ing fixed velocity pattern could re produce the ex 
cess rotation velocities found by iBeck fc Schou 
(J2000. ) . In projecting the vector velocities onto the 
line-of-sight the function sin (where (j> is the he- 
liographic longitude relative to the central merid- 
ian) multiplies the longitudinal flow velocities near 
the equator. Since the flows are largely horizontal 
the sin (j) multiplier effectively pushes the peaks in 
the Doppler pattern away from the central merid- 
ian and makes the pattern appear to rotate more 
rapidly. 

fact, 



In 



Schou 



([20031) largely removed the 
line-of-sight projection effects from the equatorial 
Doppler data by dividing the data by a function 
that approximated the function sin </> and found 
that rotation velocities were much more in line 
with those from global helioseismology (but noted 
that there were still motions relative to this that 
suggested wave-like properties for supergranules) . 
Recently, iHathawav et al. I ( 20101 ) found that the 



rotation profiles as functions of latitude deter- 
mined by the cross-correlation technique used by 
Duvall I ( 19801 ) for time lags from 2-hr to 16-hr 
could be reproduced by cellular patterns that are 
advected by a differential rotation with a peak 
velocity consistent with that found in the Sun's 
surface shear layer by global helioseismology. 

Here we measure the rotation of the pattern 
of supergranules by analy zing the same data used 
by iBeck fc SchouT (120001) and by ISchou I (l2003h . 
We execute a series of 2D Fourier transform anal- 
yses. We remove the line-of-sight projec tion ef- 
fects near the equator as was done by ISchou I 
(^2003 ) and repeat th e Fourier analysis done by 
Beck fc Schou I ( 20001) to show that the equato- 
rial rotation rate of supergranules as a function of 
wavelength matches the equatorial rotation rate 
as a function o f depth determined from global 



helioseismology (jSchou et al. 1 119981 ). This "de- 



projection" can only be done at the equator and 
is only approximate since the fiows are not purely 
horizontal. We determine the rotation rates at 
other latitudes by repeating the Fourier transform 
analysis on the raw Doppler data (without the re- 
moval of projections effects) and using a data sim- 
ulation to support our conclusions. 

2. DATA PREPARATION 

The data consist of 1024^ pixel images of the 
line-of-sight velocity determined from the Doppler 
shift of a spectral line due to the trace element 
nickel in the solar atmosphere. The images are 
acquired at a 1 min cadence and cover the full vis- 
ible disk of the Sun. We average the data over 31 
min with a Gaussian weighting function which fil- 
ters out variations on time scales less than about 
16 min, and sample that data at 15 min intervals. 
We then map these temporally filtered images onto 
a 1024^ grid in heliographic latitude from pole to 
pole and in longitude ±90° from the central merid- 
ian (Figure 1). This mapping accounts for the po- 
sition angle of the Sun's rotation axis relative to 
the imaging CCD and the tilt angle of the Sun's 
rotation axis toward or away from the spacecraft. 
Both of these angles include modification in line 
with the most recent determinations of the orien- 



tation of the Sun's rotation axis (Beck fc Giles 



2005t IHathawav fc Rightmire II2OIOI) . We analyze 
data obtained during a 60 day period of continu- 




latitudinal, 9, variation separated from the wave- 
length, A, variation such that 



n(e,\)/2TT = f{e)[i + g{X)] 



(1) 



with 



and 



/(6») = 454 - 51 sin^ 61-92 sin'' 9 nHz (2) 



g(A) = 0.045 tanh — 



2.3 — tanh 



(A -65) 



20 



/3.3 

(3) 



where the wavelength. A, is given in Mm. 



Fig. 1. — Hcliograpliic map details of the linc-of- 
sight (Doppler) velocity from SOHO/MDI (top) 
and from the data simulation (bottom). Each map 
detail extends 90° in longitude from the central 
meridian on the left and about 35° in latitude from 
the equator (the thick horizontal line). The mot- 
tled pattern is the Doppler signal (blue for blue 
shifts and red for red shifts) due to the supergran- 
ules. 

ous coverage in 1996 from May 24 to July 22. 

Line-of-sight projection effects influence the 
results so we also generate and analyze simu- 
lated data to assist in our determination of the 
actual rotation as a function of latitude and 
depth in the Sun's outer convection zone. We 
construct the simulated data from an evolving 
spectrum of vector spherical harmonics in such 
a manner as to reproduce the spatial, spectral, 
and temporal behavior of the observed cellular 
flows. The amplitudes of the spectral coefhcients 
are constraine d by matching the obse rved veloc- 



ity spectrum ( Hathawav et al. I 120001 1 with the 



radial flow component constrained by the disk 
center to limb variation in the RMS Doppler 



signal (jHathawav et al. 1 120021 ) 



The cells are given finite lifetimes and made 
to rotate by adding changes to the phases of the 
spectral coefficients (JHathawav et al. 1120101 ). The 
rotation rates are constrained by matching the ob- 
served rotation rates as functions of latitude and 
wavelength. 

The rotation rates of the cells in the simulated 
data are given by a fairly simple function with the 



3. EQUATORIAL ROTATION RATE 

We determine the equatorial rotation rate 
of the supergranules by 2D Fourier transforms 
with and without removing line-of-sight projec- 
tion effects. The line-of-sight projection effects 
ca n be m inimized using the method described 
by ISchou I ([2003) ■ The mapped Doppler veloci- 
ties near the equator are divided by a function, 
sg'n.{<j))^/{s\T? (f) + 0.01), which approximates the 
geometric factor, sin 0, that multiplies the longi- 
tudinal velocity in producing the Doppler signal. 
This "de-projected" signal and the raw Doppler 
signal are both then apodized near the limb and 
then Fourier transformed over longitude for the 
50 latitude positions that straddle the equator. 
These spectral coefficients are then Fourier trans- 
formed in time over six 10-day intervals. 

The rotation rate as a function of wavelength is 
determined by first finding the temporal frequency 
of the centroid of the spectral power for that wave- 
length. This temporal frequency is then divided by 
the longitudinal wavenumber to give the synodic 
rotation rate which is then converted to a sidereal 
rotation rate by adding a correction based on the 
rate of change of the ecliptic longitude during the 
observations (^1° day~^). 

The results of these analyses are shown in Fig. 
2 along with the rotation rate wit h depth from 
a glob al helioseismology analysis by ISchou et al. 
( 19981 ). The sidereal rotation rate of the de- 



projected supergranules as a function of longitu- 
dinal wavelength matches the rotation rate as a 
function of depth through the outer half of the 
solar convection zone and the simulation matches 
the MDI observations. The raw Doppler data give 
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Fig. 2. — The equatorial rotation rate for the 
Doppler pattern as a function of longitudinal 
wavelength. Results for the raw data are shown 
in the upper panel. Results for the de-projected 
data are shown in the lower panel. MDI results 
are shown by the black dots (with 2a error bars for 
wavelengths longer than 20 Mm). Simulated data 
results are shown by open circles. The equatorial 
rotation rate as a function of dept h from global 



helioseismology (jSchou et al. 1 119981 ) is shown by 



the large red dots while the rotation profile used 
in the simulation is shown by the solid black lines. 



faster rotation rates than the de-projected data at 
all wavelengths (in both the MDI and the simu- 
lated data) with larger increases for larger cells. 

The rotation rates seen with both de-projected 
datasets fall slightly below both the helioseismol- 
ogy results and the input profile for the simulated 
data. Experiments with the simulated data sug- 
gest that this can be attributed to the presence of 
a small radial flow component. This component is 
projected into the line-of-sight along the equator 
by multiplying by a projection factor cos (j). When 
the Doppler signal is de-projected by dividing by 
sin (j) this small signal can become large and influ- 
ence the results in this manner. 

We also see that the rotation rates for the raw 
simulated data fall systematically below the MDI 
results for wavelengths > 50 Mm. This too may 
be a result of the radial flows. 

These results do however lead to a key conclu- 



sion - that supergranules with sizes from 10 Mm 
to 100 Mm are advected by flows within the con- 
vection zone at depths equal to their widths. 

4. ROTATION PROFILES 

We determine the rotation rate of the super- 
granules as functions of latitude and wavelength 
(depth) by repeating the 2D Fourier transforms 
on the raw Doppler data for a series of latitude 
strips. Each strip is 11 pixels or ~ 2° high in 
latitude and offset from the previous strip by 5 
pixels. We repeat the procedure for 172 positions 
between ±75° latitude and for each of six 10-day 
intervals for both MDI and simulated data. Lat- 
itudinal rotation profiles are obtained for a series 
of cell wavelengths by averaging the profiles for 
all waveumbers that produce wavelengths within 
5 Mm of the target wavelength. These profiles are 
then averaged between hemispheres and smoothed 
with a 9-point binomial smoothing kernal which 
then limits the data to ±70° latitude. The results 
for averages from the MDI datasets and the from 
simulated data are shown in Fig. 3. 

The rotation velocity appears to increase at all 
latitudes with increasing wavelength. However, 
the small increase at 70 Mm wavelength (in partic- 
ular near the equator) can be attributed to projec- 
tion effects. The actual rotation rate of the con- 
vection cells in the simulation is given by equa- 
tions 1-3 and the solid lines in Fig 2. which give a 
rotation rate that increases with wavelength to a 
maximum at ~ 50 Mm at all latitudes. 

The smallest cells, 10 ± 5 Mm, are difficult to 
resolve at high latitudes but clearly show slower 
rotation than larger cells over their observed lati- 
tude range. Cells with wavelengths of 20 ± 5 Mm 
can be resolved at all latitudes and have signifi- 
cantly lower rotation rates than the larger cells. 

While this method with the raw Doppler data 
is subject to systematic offsets due to line-of-sight 
projection effects, the same offsets are present in 
the simulated data and are much smaller for the 
smaller cells. We conclude that the surface shear 
layer, in which the rotation rate increases inward, 
extends to latitudes of at least 70°. 
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Fig. 3. — The differential rotation profiles (aver- 
age longitudinal velocity relative to a frame of ref- 
erence rotating at the Carrington sidereal rate of 
456 nHz) for the cells with 30 Mm wavelengths are 
shown in the upper panel with a solid line for the 
MDI data, a dashed line for the simulated data, 
and 2tT error limits on the MDI data indicated by 
thin solid lines. The differences between the dif- 
ferential rotation at the other wavelengths from 
that found with the 30 Mm wavelength cells are 
shown in the lower panel using colors for the dif- 
ferent wavelengths - violet for 10 Mm, blue for 20 
Mm, green for 50 Mm, and red for 70 Mm. 

5. CONCLUSIONS 

We conclude that supergranules are anchored or 
steered in the subsurface flows at depths equal to 
their wavelengths. This is a simple explanation for 
the match with the rotation rate with depth from 
global helioiseismology (Figure 2 lower panel) and 
is based on well known physics. This is consis- 
tent with numerical simulations of c onvection in 
the ou termost 16 Mm of the Sun by IStein et al. 



(|201l[ ) who show that flow structures at different 
depths have diameters about equal to the depth 
itself. 

This conclusion is, however, somewhat sur- 
prising given the much smaller estimates for the 
depth of typical supergranules previously deter- 
mined from the visibility of their in t ernal flows 
using local helioseismology. lOuvall I (|l998r ) 



ules while IZhao fc Kosovichev I (l2003|) estimated 
a depth of 15 Mm. This suggests that local he- 
lioseismology is less sensitive to these deeper (and 
slower) flows and that this new method of probing 
the convection zone with supergranules can probe 
flows at greater depths. 

We also conclude that the surface shear layer 
extends to a depth of ^ 50 Mm at all latitudes. 
The increase in rotation rate with depth has long 
been suggested by observations and is attributed 
to the conservation of angular momentum for fluid 
elements movi ng inward and outward in the near 



surfa ce layers (jFoukal fc Jokipii II1975I : iHathawav 
19821 ). Measurements of this rotation rate in- 



crease from helioseisniology (jSchou et al. I Il998t 



esti- 



mated a depth of 8 Mm for typical supergran- 



Corbard fc Thompson Il2002f ) indicate that it fol- 
lows this critical gradient to depths of 10-15 Mm 
and reaches a maximum rotation rate at depths 
of 35-50 Mm. However, many helioseismology re- 
sults also suggest that the shear layer disappears 
at latitudes above about 50°. The results reported 
here, using supergranules, indicate that the shear 
layer extends to the highest latitudes probed with 
this MDI data- -70°. 

The rotation increase with depth given by 
Equation 3 follows the critical gradient (with 
dlnfl/dlnr = —1) given by angular momentum 
mixing near the surface but then drops below that 
value at greater depths as the cell turn-over times 
become longer and the convective flows adjust to 
the solar rotation. 

Global helioseismology can measure the inter- 
nal rotation rate to great depths but it gives less 
reliable results at high latitudes. Local helioseis- 
mology can measure the non-axisymmetric flows 
(as well as the axisymmetric meridional flow) but 
only in the near surface layers. Using supergran- 
ules of different sizes to probe the flows in the 
Sun's convection zone extends these measurements 
to greater depths and higher latitudes. This new 
method of probing solar convection zone dynam- 
ics should provide information complementary to 
that obtained with helioseismology. 
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